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The gluon propagator is calculated in the instanton background in a 
form appropriate for extracting the momentum dependent gluon mass. In 
background-^-gauge we get for the mass 400MeV for small p'^ independent of 
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1. Introduction 



In the following we will calculate the gluon propagator for different gauges. In 
background ^ = 1 gauge it has the simple form S"^!^ = 5"'^g^y/{jP' — M{pY). In 
section 2 we will extract the inverse propagator from the terms quadratic in the 
fluctuations around a background field. In section 3 we will average this expression 
over relevant background fields and expand it for large momentum in a way to get the 
gluon mass M{p). In section 4 we explicitly calculate Mi = M{p = oo) = 420MeV 
in the instanton background. We will also see that there are further terms which do 
not vanish at large momentum but are more difficult to calculate. To get reliable 
results for large as well as for small momentum we make in section 5 an cluster 
expansion in the instanton density. For this purpose we need the gluon propagator 
in the 1 instanton background. The relevant formulas to construct this propagator 
are listed in section 6. Because it is a bit lengthy we will restrict ourself in section 7 
to the case of small momentum. Till now as far as I know the instanton background 
is only treated at the classical level except for a calculation of the one loop action 
in the instanton background. This work should be seen as a first step to calculate 
gluonic quantum fluctuations around the instanton background and the simplest 
object to play with is the propagator. For high energies there are tools like operator 
product expansion with instantons included indirectly in the condensates 0. This 
is another reason for the restriction to small momentum. 

2. Gluon Propagator 

The first task to obtain a formal expression for the gluon propagator in a background 
field is to expand Cqcd[-A + B] in the fluctuations B'^ around our background A^. 
The term quadratic in 5° is then by definition the inverse gluon propagator. For 
A"^ we will later use our instanton gas. We will use the QCD-Lagrangian 



where we have rescaled the fields in such a way that the coupling-constant- 
dependence is in front of C. We will entirely work in Euclidian space with metric 
Sfj^u instead of gf^i, because instantons do not make sense in Minkowski space. At the 
very end we can simply rotate back to Minkowski space. With these conventions 
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where D'^ = d^Sab + facb^^^ is the covariant derivative with A inserted instead of A; 
similarly G^j, = As usual the terms which are total derivatives disappear 

after integrating C We see that if A solves the equation of motion the linear term 
vanishes, as it should be. If we choose background gauge DlfB'^ = the last term 
of the 0(i?^)-contribution vanishes. Now we can read the inverse gluon propagator 
from the terms quadratic in B (from now on we will omit the bars over A, G and 
D because the unbared objects won't be needed further): 

= ^{-D''pDf5^u - 2facbGl^) . (2) 

We will also omit the l/g"^ in front of the propagator which is a result of the rescaling 
of fields anyway. For further manipulations some abbreviations are useful: 

{F'')ab = tfacb , [F\F']=ifabcF' , tr.F^F' = NJ^" , N, = 3 , 
{P,t = iDf , P, = id, , P, = p, + A, , 

P^X=[p^,X]+Xp^ = i{d^X)+Xp^ . (3) 

The last equation has only been quoted to show that p and P will be used in 
operator sense. F"- are the generators in adjoint representation and fabc are the 
structure constants of the color gauge group SU{Nc). With these abbreviations we 
can now write 

S% = 5,u/P'' , V,, = {A'+p-A + A-p)6,, + 2tG^, . (4) 

Sq is the free gluon propator with no background and V can be interpreted as 
an interaction potential caused by the background. The QCD-Lagrangian in the 
background |I| can be written in another form more suitable for ordinary perturbation 
theory: 

g'CQCD = + Ib^V^^^B^ + UcfaedB^BtB^^A^^ (5) 

The terms independent and linear in B has been omitted because for a given back- 
ground the first term is irrelevant and the second can be deleted by shifting B. To 
determine the weight of a specific background they are of course the most important 
terms. The Feynman graphs of the additional terms (with g recovered) are depicted 
in figure . It should be noted that for small coupling constant g the second graph 
can be treated perturbatively but not the first. So our main concentration lies on 
the first term. 
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Figure 1: Feynman rules for gluons in a background 



3. Propagator in Statistical Background 

The next step is to use some approximation scheme to calculate the propagator 

S={S^^ + V)-^ = S^{\ + T)-^ , T = VSo . (6) 

For large momentum p, 5*0 and therefore T are small and we can expand S in powers 
of T: 

S = Sq{\-T + - T^ + ...) . (7) 

Note that S{x,y) = {x\S\y) is generally not invariant under translations and ro- 
tations because the background and therefore T are not. Actually we are not 
interested in the propagator for a particular background configuration, but only in 
the average over all relevant configurations. Here we do not mean the functional 
integration over quantum fluctuations around the empty vacuum, but fields other 
than the perturbative ^4^ = -vacuum which minimize the action / Cdx. 

S = So{l-T + T^-T^+...) , (8) 

where the bar denotes averaging over relevant configurations. Details are specified 
below. If the background is statistically invariant under translations then S is trans- 
lationally invariant and therefore diagonal in momentum space S{p,q) = {p\S\q) = 
S{p)^{p — q). If we would evaluate the terms in the series we would get an expansion 
of S in the form 

Sip) = l/p' + ci/p^ + C2// + . . . (9) 

and the pole remains at zero — and gets even worse with higher terms. What we 
want is S" in a form like S{p) = [p^ + M{pY)~^ with M{p) bounded for large and 
small p and interpreted as momentum dependend gluon mass0. So let's invert (§) 

= p2 + M{pf = (1-T + T2-T3 + .. .)-^Sq^ (10) 

^ The " + " in front of M will change to the more familiar " — " when we rotate back from 
Euclidian space to Minkowski space. 
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and expand it once again in T. Without averaging this would just be a geometrical 
series which, expanded, would give the original formula S^^ = Sq^ + V. With 
averaging the different terms are now in no relation and expanding and sorting with 
respect to powers of T yields 

S'^ = (t + T-{T^-T^) + {T3-TT^-T^T + T^))S^^ + 0{T^) (11) 
= So^ + V-{VSoV-VSoV) + ... 

= Sq^ + M\ M"^ = Ml- Ml + Ml = V, MI = VSqV -VSqV. 

In the next section we introduce the concepts of instanton gas calculation to deter- 
mine Ml, which is actually very simple. 



4. First Order Gluon Mass 

The classical approximation to a quantum theory in the language of path integrals 
is to consider only the configurations which minimize the classical Euclidian action 
5*. In QED there is only one minimum, namely = 0, but in QCD there are other 
local minima, called instantons: 

A%{x) = Ofvi. _ I'^^l _ I'^y ^ ^2 ' (12) 

V^^^u = Va^lu = - ^eabcebcf^u for instantons , 

vi^^u = "la^iu = ea^ui + ]j^(^abc(-bc^lu for anti-iustantous . 

A'j^ is an (anti)instanton of size p at position zj in singular gauge and orientation 
O/, where Oj is a rotation matrix in the adjoint representation of SU{Nc). rja^u are 
the 't Hooft-symbols with some properties given in the appendix of 0. Whereas 
the single instanton is an exact solution of the equation of motion D^^G''^^ = 0, the 
multi-instanton- field A = J2i Aj minimizes S only approximatly, but nevertheless 
gives a significant contribution to the functional integral. Quantum fluctuations 
around the instanton renormalize the action 5* = Svr^/^f^ and make small instantons 
less important because the QCD coupling constant g decreases for small distances 
0. Repulsion between instantons do not allow too large instantons 0, 0, A]. So 
there is a narrow region of allowed values for the instanton radius and I will use an 
average radius p in my calculations. The second parameter we need is the instanton 
density n which must be determined experimentally from the gluon condensate . 
The ratio Lq/p is calculated in p, |Q: 

n = iV/\/4 = l/I^^ = ^(G^.Gn = (200MeV)L,. , Vp = 3.2,,eo.. (13) 
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Now we plug the N-(anti)instanton-field 



N 



A = J2Aj (14) 



1=1 



in our expressions for the gluon propagator ([TT|) . To do this we must average some 
functions f{A) of A like in V. We will use the instanton gas approximation without 
interactions between instantons. The whole effects of interactions are summarized 
in the values of Lq and p. So the position and orientation are equally distributed 
and independent for different instantons: 



^ 1 /• r 

f{A) = J[-Jd%JdOjfiA) , (15) 



where J dzj averages over the position of instanton I contained in a large box of 
volume V4. / dOj is the Haar- measure and averages over the group of orientations 



d01 = l , 0^^ = J dOO^^ = , 

O^bQ^ = fdO O'^^O^'^ = — (16) 
J N^-l 

We now can easily calculate V^i, = (A^ + p ■ A + A ■ p)5^y + 2iG^y : 



A^^ = NAf^ = NO^b ■ . . . = and G"^^ = (17) 

because is antisymmetric in the Lorentz indices and Lorentz invariant^ after 
averaging, which is impossible. 

N 

A^ = Y.Aj + Y.A;A; = NAj + N{N-l)AjAj = NAj , 
1=1 ij^j 

where we have used the fact that different instantons are independent and (|T7|). 
Inserting now the specific form of the instanton we yield 



{A^)ab = y^JdZlJ dOifacdA%lfdebA% 
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-j^ -^np dab , ^ ~ y' ~ instanton density (18) 
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^Precisly we should say 0(4) invariant because we are in Euclidian space. 
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where we have used facdfdcb = -NcSab and ■q^^^'ql^^p = Z5„p. Inserting this into 
Mf = F = we get 



127r2/V 



Ml = ^ ^^^^(p/Lo)i^o ~ = 420 MeV (19) 

Unfortunately M2 also contains terms which survive the large p limit. To see this 

one can count the number of "p's occuring in M| which will give us the dominant 
behaviour of M| for large p. VSqV = Mf/p^ vanishes for large p but in principle 
V contains a p in the nominater {pA) and Sq = l/p"^ and V SqV could be finite for 
large p. To be more definite consider 

VSqV = ipA + Ap)SQ{r)A + Ap) + other terms = AA^^^J^A^ + . . . 

where we have used [p, Ai] = id^Aj — 0. 

{x\VSoV\y) = 4A,{x){x\^\y)A,{y) + ... 

Prom 5^^{x\p^p^/p^\y) = {x\y) = 5{x - y)we can conclude that(x|p^p^/p^|y) = 
\5fj,i,S{x — y) + traceless terms. 

{x\VSoV\y) = A^{x)A^{x)6{x - y) + . . . = {x\A'\y) + ... 

So M| = A"^ + . . . = Ml + . . . contains a term which fully cancels calculated 
above. In the next section we use an approximation which overcomes this problem 
but actually restrict ourself in a calculation for small p. 



5. Propagator in statistical Background (Expansion in Instanton Density) 

In this section we will make an expansion of the gluon propagator in the instanton 
density n which will be valid for all euclidian momenta p especially for small p. 

The average of an arbitrary power of the 1-instanton field is proportional to the 
inverse volume AJ y for n > 1. So, for example, the expansion of the square 
of the multi-instanton configuration in the instanton density is 



N 




and more general 

= NA^ + 0(n^) , W = NVp + 0(n^) for n > 1, (21) 
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where Vj = V{Aj) defined in (|) with A replaced by Aj and similar for T. Let us 



now sum up all terms in (|TT]) linear in n: 



S ^ = {l + T-T^ + T^-...)So^ + 0{n^) 



[1 + N{Tj - T] + Tf-.. .))S^^ + 0{n 



_ _^ 

= {\ + NT,jf)S,' + 0{n') = S,' + NV,ff + 0{n^) (22) 

Te// = Ti-Tf + Tf-... = Tj-Ti{Tj-Tf + ...)= Tj-TiT,ff 

Veff = Vl~VlSoVeff =^ Veff = S^\So~Sl)S^' with (23) 

= Sq^ + Vj is the propagator in the 1-instanton background. (24) 

More generally we can make a cluster expansion of an arbitrary function of A in the 
following form 

N 

f{A, + ... + A^) = f{A) = ^(r) E i-y'%)fiAi + ... + A,)= (25) 

/=0 k=0..l 



= /(O) + NfiA,) - /(O) +- N{N - l)f{A, + A2) - 2f{A,) + /(O) +0{n') 

, ■■ L , ' 

0(1) 0(n) 0(n2) 

where the first line is an identity even without averaging. Inserting a Taylor expan- 
sion for / in the second line and using the indistinguishability of different instantons 
one can see that all monomials in the fc-th term contain more or equal than h differ- 
ent instantons. So the average will factorize in k factors each proportional to n and 
therefore the fc-th term is indeed proportional to ^ . It is easy to generalize (25) for 
two or more species of fields. If A is a field of Ni instantons and Nj antiinstantons 
we get to first order in n 



f{A) = /(O) + NjfiAj) - /(O) + NjfiAr) - /(O) + 0{n'). (26) 
If we insert the propagator S in (|2^) we get 



S = S{A) = S{0) + NS{Ai) - S{0) + 0{n^) = So + NSi - ^0+ 0{n^) (27) 
which is after inversion up to O(n^) just (^2]) . 



6. QCD Propagators 

In section 5 we have seen that it is enough to know the 1-instanton propagator to 
calculate S to first order in n. Luckily this propagator is known, unfortunatly it is 
a bit lengthy and suffers from a divergency. 

The gluon propag ator Sf^^ with spin 5* = 1 in adjoint colorQ representation (C = 1) 
can be constructed out of the ghost propagator A"'' (5* = 0, C = 1) which is explicitly 

^ Often called isospin Q 
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known in the 1-instanton background. The general formulas how to construct a 
propagator of given spin S out of the corresponding scalar propagator with same 
color C in a selfdual background are shown below. They are derived and more 
thoroughly discussed in f^. 

Notations: 

= T'^Al , [T^T^] = ze„,eT^ , = generator of SU{2), 

r in scalar (C = 0) \ 

T'^ = < r"/2 in fundamental {C = \) > representation (28) 
[ ie.a- in adjoint (C = 1) J 

Spin propagator A : 

A-i = p2 ^ ^ _ p-2 (29) 
I propagator S: (not used hut only stated for completeness) 
S-' = r=^^P'^ =^ 5 = ^Ai^ + Afi^ for G^, = G^, (30) 



Spin 1 Propagator S: 

S-^ = P'^K,, + 2zG^, - (1 - i)P^P, ^ 

5^. = gM-p^PpA^P, - (1 - OP^A^P, for G^, = G^,, (31) 

There are some comments in order. The above formulas are valid for an arbitrary 
color representation, we will need them only for C = 1. For S ^ 1 there are zero- 
modes and the propagator is only the inverse of the kernel in a subspace orthogonal 
to the zeromodes. The implications and problems will be discussed in section 8 
when they show up explicitly. The Spin 1 kernel is the quadratic term of the QCD- 
Langrangian (|l]) with gauge fixing term ^{D^^^'B^^Y in slight generalisation to the 
^ = 1 case considered previously. 

With 

n(x) = l + 4 , F{x,y) = l + p'^^^^ , (32) 
Tp = (r, i) , tI = (r, -i) , r^rj = 6^^ + i^puTa 
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we can write the ghost propagator for 1 instanton in the form 



,a6/^ .A ^ |tr TaF{x,y)nF{y,x) _ 



47r2(x - yyU{x)U{y) 

_ Sab p'^Sgb I '^p'^egbcVc^.uXf.yi, 

ATi'^{x — yY 47r2(x2 + p2)(?/2 + p2) A.tt'^^x — yY^x"^ + p'^){ii'^ + p^) 
^ 2p'^{{{xyf - x'^y'^)Sab + egbc'gci.vXf.yy + Vq^uVbpqXf.y^Xpy^ 
47r2(x — y)2(a;2 + p2)a;2(l/2 + p2)?/2 

For simphcity we have placed the instanton at the origin zi = in standard orien- 
tation O"'' = 5"^. It is possible to write down the gluon propagator Sj'^^, explicitly 
in which we are finally interested in using (^TJ) but the expression will be rather 
lengthy. While Aj'' can be averaged and represented in momentum space exactly 
with the help of modified Besselfunctions this seems not to be possible for Sj'^^^ and 
we have to expand 5* for large and/or small momenta (or work much harder to solve 
the complicated integrals in terms of special functions). So we will never use the 
full expression for S. 



7. Propagators for small momentum 

The calculation simplifies significantly if we expand A(p) or S{p) for small momen- 
tum. Because p always occurs in the dimensionless quantity (pp) the lowest order 
in p can be found by keeping only terms of lowest order in p or differently stated: 
Small p corresponds to large x and for x 3> p p is neglegible. To warm up let us 
start with A from (P3|): 

Af = A^^ - p^W^ + 0{p^) (34) 



Atx'^[x — yY 

J ah I \ Sah '^^ahc^cji.uX ^y^ 



47r2x2y2 ^t:'^{x — yY'X^y'^ 

Ao(p) = is the free (A^ = 0) ghost propagator. After reintroducing the in- 
stanton position z we can now average Aq — A/. The e-term in W will be killed by 
averaging over the instanton orientation: 

1 ( dh 



(.|A?-Af|,) = 4.W^/j;5(737)3S3(73i;)i . (35) 

The last integral is infrared diverent but noticing that Xj^-n'^ix — yY is the fourier 
transformation of l/p^ we can rewrite the above expression in the form 

VM'^^i\y) = 4vrV / d^z{x\\\z){z\-^\y) = 47rV(x|^|i/) 

J p p p 
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The divergence is now hidden in the fact that the coordinate representation of 
does not exist. In fact we are not interested in the coordinate representation but 
in the momentum representation and so the divergence is spurious. Inserting A in 
(^) we get 

Vf;; = Ao ^ Ao ^ = P'^^y = y^P' = (36) 

and finally 

= 0) = We// = An'p'n ^ (420MeV)2 
which has to be inserted in the ghost propagator Aab = Sabip^ + ^ghostip))^^ ■ 

What can we learn from this? At first M ^ 0. A scalar particle moving in the 
instanton background gets a dynamical soft mass. The pole at p = has vanished. 
Secondly M{p 0) is finite, reasonable to think that M is bounded (has no poles) 
at least in the Euclidian region p"^ > 0. From an exact averaging of (^) one can 
see that M is finit for all momenta This was our motivation to prefer the 

representation (|Tl|) instead of (P). Let me make a last comment on this inversion. 
In ordinary perturbation theory one calculates the one loop self energy, sums up 
the series of one loop two consecutive loops and so on to bring the selfenergy in the 
denominator of the propagator which actually is equivalent to a simple inversion 
and hopes to improve the result with this partial resummation. The arguments for 
doing this are not better or weaker than here - they are essentially the same. 

Let us now calculate the gluon mass at zero momentum along the same lines. To 
do this we must insert (0) into (pT)). Expanding also up to order 

P,=P, + p'A, + Oip') , A, = ^^"^7-^- 

we get 

P^A'jP, = p^Aop, - p%iAoW + WAo)p, + p^p^AoA, + A^Aop,) + O(p^) 

The free gluon propagator in i?^-gauge is well known or can be obtained from (|3lD 
setting Gfj,i, = 0: 

= %upaPp^lPa - (1 - OPp^IPv = ^{^t^u - (1 - 0^^) 

S'u - = p'[q,.paPpi^oW + WAo)p^ - (1 - OPm(AoW^ + WA^ 



5°. - = 2p'w{p%, - (1 - Op,p^) = ^pVwsl 

where we have used in the last line that p comutes with averaged quantities like W. 
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Maybe the relation Mgi^on ~ '^^ghost is valid even for larger p. Till now we have 
dealt with 2-color QCD = 2. The extension of our result to the real world needs 
a thorough adaption of all formulas in this work and P] to A^^^ = 3. After averaging 
there will appear several factors Nc/{N^ — 1) like in ( |T9| ) and 3/{N^ — 1) compared 
to 2/3 and 1 in the A^^^ = 2-case. So most probable M^^^^^ has to be multiplied with 
a factor between 3/8 and 1/2 yielding a gluon mass of approximately 400MeV for 
A^. = 3. 



8. Conclusions and further developments 

In our whole calculation we have ignored the zeromodes. All gluon field fluctua- 
tions must be orthogonal to these which is achieved by using a gluon propagator 
orthogonal to the zeromodes or otherwise stated by subtracting out from (|3T|) the 
projections on the zeromode subspace. This procedure also deletes a divergence 
comming from the second term in ( p3D squared. There are some further problems 
to be solved because the scalar product of the zeromodes with the propagator does 
not exist [^. But all this concerns a term which is proportional to which was 
irrelevant in our zero momentum approximation. Also the principle mechanism of 
mass generation does not depend on the zeromodes which can be seen from the 
ghost propagator because in the spin case there are no zeromodes. This should be 
contrasted with the fermionic case where the zeromode are the main ingredient for 
mass generation and the so called zeromode approximation simplifies calculations 
enormously. 

The next step may be the calculation of some gauge invariant correlation function 
like a glueball correlator or topological susceptibility. As usual this involves products 
and integrals over propagators but with the difference that we have to use ( PTD 
instead of the simple free propagator. Without further simplification this may cause 
a headache. 
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